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Method of improving fingerprint images 



The invention relates to a method of improving fingerprint images, wherein 
recorded images are processed on image areas of the fingerprint images by using Gabor 
filters. 

Of the methods currently used to allow clear identification of different 
5 fingerprints with maximum reliability, most work on the basis of extraction of minutiae and 

comparison of position and, where applicable, nature and direction of the minutiae found. 

Minutiae are branches and ends of individual grooves in the fingerprint. The minutiae are 

extracted from a given fingerprint using known image processing methods. Errors in the 

production of fingerprint images hinder the detection of rninutiae. Such errors can, for 
10 example, arise from different composition of the finger surface at the time of photographing, 

for example, sweat, dirt, dryness and the like. 

US 5.659.626 discloses a fingerprint identification system in which processing 

is performed in accordance with the type of the main claim using Gabor filters. Gabor filters 

are described further in Lin Hong, Yifei Wan, and Anil Jain: "Fingerprint Image 
15 Enhancement: Algorithm and Performance Evaluation", IEEE Transactions on Pattern 

Analysis and Machine Intelligence, Vol. 20. 8, August 1998. 

These methods indeed give an essential improvement in fingerprint images for 

subsequent image analysis. It has, however, been found that the results are not always good 

for all variants of fingerprints and parts thereof. 
20 Therefore, according to the invention, to further improve fingerprint images, 

the Gabor filters are adapted to characteristic quantities of the image in the relevant image 

area to be processed. 

In order to be able to extract suitable characteristic quantities from the content 

of the image areas for adaptation of parameters, the method according to the invention is 
25 preferably characterized in that the image areas are square and have a side length at which the 

spacing and direction changes of the grooves are relatively small but comprise several 

grooves, such that statistical mean values and the standard deviations can be formed for the 

spacings and directions. 
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An essential improvement in fingerprint images is achieved by an 
advantageous embodiment of the method according to the invention in the adaptation of the 
image area to be processed is determined by the following characteristic quantities: 
mean direction of the grooves in the image area, 
5 mean spacing of the grooves in the image area, 

change of direction of the grooves in the image area, and 
change of spacing of the grooves in the image area. 

The method according to the invention is particularly characterized in that, of 
an ellipse which encloses a Gaussian bell curve representing the Fourier transform of the 
10 Gabor filter, the length of the main axis directed towards the origin increases as the spacing 
of the grooves decreases, the length of the other main axis increases as the direction change 
increases, and the intersection of the main axes is determined by the mean direction and mean 
frequency of the grooves in the relevant image area. 

When combining the image areas distorting artifacts which are due to the 
1 5 filtration varying from image area to area may occur in the peripheral areas of the individual 
Gabor-filtered image areas, which artifacts affect further image processing. Therefore, in a 
refinement of the method according to the invention, the image areas overlap and before 
being combined, the overlapping image areas are weighted with a function decreasing 
towards the edges. It is thereby achieved that these artifacts disappear or are minimized and 
20 hence have no negative effects on the required quality of the overall Gabor-filtered image 
recombined from the individual tiles. 

These and other aspects of the invention are apparent from and will be 
elucidated with reference to the embodiments described hereinafter. 



25 

In the drawings: 

Fig. 1 shows fingerprint images before and after processing by means of the 
method according to the invention with the result of an improved groove contour, 

Fig. 2 shows fingerprint images before and after processing by means of the 
30 method according to the invention with the result of an improved groove separation, 

Fig. 3 shows a flow chart to explain the method according to the invention, 

Fig. 4 shows views to explain the combining of tiles, and 

Fig. 5 shows a plan view of the FFT-transform of a square quadratic image 

area. 
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Because of the specifications for patent drawings which state that the figures 
must not contain gray shades, the fingerprint images in Figures 1 and 2 are shown largely as 
pure black/white images. 

Fig. la shows an unprocessed fingerprint with an error point caused, for 
example, by a scratch on the surface of the scanner. The interrupted grooves are repaired 
(Fig. lb) with Gabor filters adapted optimally according to the method of the invention.. 

Fig. 2 demonstrates a second property of the method according to the 
invention: separation of the distorting wash-out of the grooves due to harmful environmental 
influences in the marked areas of the sensor image (Fig. 2a). Here, Fig. 2b shows the 
fingerprint after processing with the Gabor filters optimally adapted by means of the 
invention. 

To perform the method according to the invention in accordance with the flow 
chart shown in Fig. 3, a corresponding computer program can be run on a computer. After a 
start at 1 , the calculation of local gradients G is performed on a suitable grid of, for example, 
8x8 pixels. In a step 3, local groove spacings R are also calculated on the suitable grid. At 4, 
the fingerprint image is divided into tiles K of suitable size. A suitable size is, for example, 
32 x 32 pixels for a sensor-induced image resolution of 500 dpi. 

In program step 5, a counter c is initialized, i.e. set to zero. The number of tiles 
is also initialized with N = #K. After a query 6 whether c = N, if not applicable c is 
incremented at 7. Then, at 8, parameters are calculated for construction of the Gabor filter 
according to the characteristic quantities. Here, mean is the mean value and std is the 
standard deviation. At 9, the Gabor filter GF for the tile K(c) is constructed as a function of 
the parameters calculated in 8. 

In program step 10, the Gabor filter GF concerned is applied to the tile K(c), 
where for filtration the fast-Fourier-transformed (FFT) tile K(c) is multiplied elementary 
element with the Gabor filter GF. The result is inversely transformed (IFFT) and the real part 
forms the improved content K(c)' of tile K(c). At 1 1, the results of program step 10 are 
evaluated with a function W which will be explained in more detail below with reference to 
Fig. 4. In program step 12, the filtered and evaluated content K u (c) of tile K(c) is then 
inserted in the resulting improved fingerprint image. 
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After this program step, it is checked at 6 again whether c has reached value 
N. If this is not yet the case, method steps 7 to 12 are performed for the following tile. After 
all tiles have been processed, the program ends at 13. 

When combining the image areas, distorting artifacts, which are due to the 
filtration which differs from image area to image area, may occur in the peripheral areas of 
the individual Gabor-filtered image areas, which artifacts affect the further image processing. 
Therefore, according to Fig. 4, the image areas overlap and, before being combined, the 
overlapping image areas are weighted with a function decreasing towards the edges. It is 
thereby achieved that these artifacts disappear or are minimized and hence have no negative 
effects on the required quality of the overall Gabor-filtered image recombined from the 
individual tiles. Fig. 4a shows a tile K(c), and in dotted lines, the adjacent tiles which each 
overlap. Fig. 4b shows an example of a weighting function W as a function of distance x. 
Whereas in the central area of the tiles a maximum weighting is applied, this diminishes in 
the overlap areas up to the edges. 

In the frequency range, a Gabor filter consists of a two-dimensional Gaussian 
bell curve which, in the plan view in Fig. 5, is surrounded by an ellipse. A main axis of the 
ellipse in the case shown is in principle aligned to the center point of the filter mask. This 
main axis in the view of the filter lies in the spatial area of the normals to the wavefront. 
Thus, the four parameters of direction, direction change, frequency and frequency change 
determine the position and form of the ellipse. 

The center point of the ellipse describes the maximum of the bell curve and 
lies precisely at the intersection of frequency and direction. The frequency space is viewed as 
follows: The origin lies in the center of the Gabor filter and can be described with the 
amplitude of the equal part. Concentric circles about the origin (0,0) contain the same 
frequency in different directions, where the frequency is limited at the outer edge by the FFT- 
specific limit frequency. In this simplified view, the imaginary part of both the transform and 
the inverse transform is disregarded. Thus, the frequency space can be regarded as the polar 
co-ordinate representation of frequency and direction. 

A possible type of calculating a Gabor filter adapted according to the 
invention is explained below. The implementation has size 32x32 for the square FFT/IFFT 
tiles K. The Gabor filter in the frequency area is formed pixel by pixel for: 
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(x,y)< 



1 2 3 , (32-1) 

; 1 '- 1+ ^- 1+ i6-- 1+ ^-'- 1+ V 

GF{x,y) = txp[~v T A-^ [1] 



and as follows: 



with v 



-0- 



J_(°° smG ^ w here m# is the mean groove spacing and mG is the mean 
mR^sinmG) 

direction of the groove to the inner part of the tile being processed, calculated, for example, 
5 by means of 

mR = T77 Z P] 

where M is the quantity of elements of R existing on the tile. 

For tile size 32 equation [2] becomes 

+ 12, y + 12) + J?(* + 12, y 4- 20) + -I- 20, y + 12) + gpc + 20, y + 20) ^-j 

4 

10 with (*,;;) as co-ordinates of the top left-hand corner of the tile processed in each case. 

The value mG can be taken e.g. directly from a suitably smoothed grid of 
directions, as the innermost grid point of G in the tile being processed. This can be performed 
instead of the mean value formation and is possible because of the previously performed 
smoothing of the direction grid and prevents a complex mean value formation of the 

15 directions. 

Matrix A in [1] is derived from 



fcosmG -smmGYa? 0 V cosmG sinwG 
[sinwG cosmG \ 0 cr; 2 ){- sin mG cosroG 



3 



[4] 



where parameters a x and a y ofR, G, mR and mG can be derived e.g. as follows: 



cr = max 



A> and u y =max 



2^3-«J 



[5] 
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1 . A method of improving fingerprint images, wherein recorded images are 
processed on image areas of the fingerprint images by using Gabor filters, characterized in 
that the Gabor filters are adapted to characteristic quantities of the image in the relevant 
image area to be processed. 

5 

2. A method as claimed in claim 1 , characterized in that the image areas are 
square and have a side length at which the spacing and direction changes of the grooves are 
relatively small but comprise several grooves, such that statistical mean values and the 
standard deviations can be formed for the spacings and directions. 

10 

3. A method as claimed in any one of claim 1 or 2, characterized in that the 
adaptation of the image area to be processed is determined by the following characteristic 
quantities: 

mean direction of the grooves in the image area, 
1 5 mean spacing of the grooves in the image area, 

change of direction of the grooves in the image area, and 
change of spacing of the grooves in the image area. 

4. A method as claimed in claim 3, characterized in that, of an ellipse which 
20 encloses a Gaussian bell curve representing the Fourier transform of the Gabor filter, the 

length of the main axis directed towards the origin increases as the spacing of the grooves 
decreases, the length of the other main axis increases as the direction change increases, and 
the intersection of the main axes is determined by the mean direction and mean frequency of 
the grooves in the relevant image area. 

25 

5. A method as claimed in any one of the preceding claims, characterized in that 
the image areas overlap and before being combined, the overlapping image areas are 
weighted with a function (W) decreasing towards the edges. 
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